Maximized Orbital and Spin Kondo effects in a single-electron transistor 
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We investigate the charge fluctuations of a single-electron box (metallic grain) coupled to a lead 
via a smaller quantum dot in the Kondo regime. The most interesting aspect of this problem resides 
in the interplay between spin Kondo physics stemming from the screening of the spin of the small 
dot and orbital Kondo physics emerging when charging states of the grain with (charge) Q = 
and Q = e are almost degenerate. Combining Wilson's numerical renormalization-group method 
with perturbative scaling approaches we push forward our previous work [K. Le Hur and P. Simon, 
Phys. Rev. B 67, 201308R (2003)]. We emphasize that for symmetric and slightly asymmetric 
barriers, the strong entanglement of charge and spin flip events in this setup inevitably results in a 
non trivial stable SU(4) Kondo fixed point near the degeneracy points of the grain. By analogy with 
a small dot sandwiched between two leads, the ground state is Fermi-liquid like which considerably 
smears out the Coulomb staircase behavior and hampers the Matveev logarithmic singularity to 
arise. Most notably, the associated Kondo temperature T^"' 4 ' might be raised compared to that in 
the conductance experiments through a small quantum dot (~ IK) which makes the observation of 
our predictions a priori accessible. We discuss the robustness of the SU(4) correlated state against 
the inclusion of an external magnetic field, a deviation from the degeneracy points, particle-hole 
symmetry in the small dot, asymmetric tunnel junctions and comment on the different crossovers. 

PACS numbers: 75.20.Hr,71.27.+a,73.23.Hk 



I. INTRODUCTION 

Recently, quantum dots have attracted a considerable 
interest due to their potential applicability as single elec- 
tron transistors or as basic building blocks (qubits) in the 
fabrication of quantum computers^ In the last years, a 
great amount of work has also been devoted to studying 
the Kondo effect in mesoscopic structures^ A motivation 
for these efforts was the recent experimental observation 
of the Kondo effect in tunneling through a small quan- 
tum dot in the Kondo regime^^ In these experiments, 
the excess electronic spin of the dot acts as a magnetic 
impurity. Let us also mention that the manipulation 
of magnetic cobalt atoms on a copper surface and more 
specifically the observation of the associated Kondo res- 
onance via spectroscopy tunneling measurements^ also 
stands for a remarkable opportunity to probe spin Kondo 
physics at the mesoscopic scale but in another realm (not 
with artificial structures). 

A different set of problems relating the Kondo effect 
to the physics of quantum dots is encountered when in- 
vestigating the charge fluctuations of a large Coulomb- 
blockaded quantum dot (metallic grain) & More precisely, 
one of the most important features of a quantum dot 
is the Coulomb blockade phenomenon, i.e., as a result 
of the strong repulsion between electrons, the charge 
of a quantum dot is quantized in units of the elemen- 
tary charge e. Even a metallic dot at a micronmetric 
scale can still behave as a good single-electron transis- 
tor. When the gate voltage V g is increased the charge 



of the grain changes in a step-like manner. This be- 
havior is referred to as a Coulomb staircase. Moreover 
when the metallic dot is weakly-coupled to a bulk lead, 
so that electrons can hop from the lead to the dot and 
back, the dot charge remains to a large extent quantized. 
This quantization has been investigated with thorough- 
ness both theoretically 9, 10 ' 1 ii2i and experimentallyii 3 . It 
is important to bear in mind that such a problem is 
intrinsically connected to an orbital or charge Kondo 
effect. 9 Indeed, near the degeneracy points of the average 
charge in the grain one can effectively map the problem 
of charge fluctuations onto a (planar) two-channel Kondo 
HamiltonianiiiSiiii with the two charge configurations in 
the box playing the role of the impurity spin2*i£ and the 
physical spin of the conduction electrons acting as a pas- 
sive channel index. (This mapping is a priori valid only 
for weak tunneling junctions between the grain and the 
lead). For accessible temperatures - in general, larger 
than the level spacing of the grain - spin Kondo physics 
is not relevant^ The quantity of interest is the average 
dot charge as a function of the voltage applied to a back- 
gate. Note that the average dot charge can be measured 
with sensitivity well below a single charged Unfortu- 
nately, only some fingerprints of the two-channel Kondo 
effect were recently observed for a setting in semiconduc- 
tor quantum dots: 29 Indeed, the non-Fermi liquid nature 
of the two-channel Kondo effect is hardly accessible in 
the Matveev's setup built on semiconducting devicesi^i 
On the one hand, the charging energy of the grain must 
be large enough to maximize the Kondo temperature 
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Tk on the other hand the level spacing must be small 
enough compared to Tk- It is difficult to satisfy these 
two conflicting limits. A better chance for observing the 
two-channel Kondo behavior may be reached if tunnel- 
ing between the lead and the grain involves a resonant 
level since it offers the possibility to actually enhance the 
Kondo temperature of the system^ 




When Q=0 and Q=e are degenerate 
orbital pseudo-spin T=l/2 



FIG. 1: Schematic view of the setup. A micron scale grain (or 
large dot) is weakly coupled to a bulk lead via a small dot in 
the Kondo regime which acts as an S=l/2 spin impurity. The 
charges of the grain and the small dot are controlled by the 
gate voltage V a and Vd respectively. The auxiliary voltages 
can be used to adjust the tunnel junctions. 

In this paper, the setup we analyze consists of a single- 
electron box or grain coupled to a reservoir through a 
smaller dot (Figure 1). We assume that the smaller 
dot contains an odd number of electrons and eventu- 
ally acts as an S=l/2 Kondo impurity. 2 Typically, when 
only charge Kondo flips are involved the low energy 
physics near the degeneracy points is well described by 
a two-channel Kondo model, in particular the capaci- 
tance peaks of the grain exhibit at zero temperature a 
logarithmic singularity at the degeneracy points which 
ensures a nice Coulomb staircase even for not too weak 
couplings between the quantum box and the lead£ In 
our setup, the Kondo effect has now two possible origins: 
the spin due to the presence of the small dot playing 
the role of an S=l/2 spin impurity and the orbital de- 
generacy on the grain. Combining Wilson's numerical 
renormalization-group (NRG) method with perturbative 
scaling approaches we extend our previous workjS 3 - and 
emphasize that at (and near) the degeneracy points of 
the grain, the two Kondo effects can be intertwined. The 
orbital degrees of freedom of the grain become strongly 
entangled with the spin degrees of freedom of the small 
dot resulting in a stable fixed point with an SU(4) sym- 
metry. This requires symmetric or slightly asymmetric 
tunneling junctions. Furthermore, the low energy fixed 
point is a Fermi- liquid which considerably smears out the 
Coulomb staircase behavior and hampers the Matveev 
logarithmic singularity to arise^ Remember that the ma- 



jor consequence of this enlarged symmetry in our setup 
is that the grain capacitance exhibits instead of a loga- 
rithmic singularity, a strongly reduced peak as a func- 
tion of the back-gate voltage, smearing charging effects 
in the grain considerably. It is also worth noting that the 
Kondo effect is maximized when both Kondo effects oc- 
cur simultaneously. In particular, the associated Kondo 
temperature T K U ^ can be strongly enhanced compared 
to that of the Matveev's original setup which may guar- 
antee the verification of our predictions. We stress that 
the Coulomb staircase behavior becomes smeared out al- 
ready in the weak tunneling limit due to the appearance 
of spin-flip assisted tunneling. A different limit where 
the small dot rather acts as a resonant level close to the 
Fermi level has been studied in Refs. 22 24] where in 
contrast it was shown that the resonant level has only a 
slight influence on the smearing of the Coulomb blockade 
even if the transmission coefficient through the impurity 
is one at resonance. This differs markedly from the case 
of an energy-independent transmission coefficient where 
the Coulomb staircase is completely destroyed for perfect 
transmission^^ Furthermore, the charge of the grain in 
such a device can be used to measure the occupation of 
the dot The resonant-level behavior of Ref. [24J is also 
recovered in our setup when an orbital magnetic field is 
applied. 

Let us mention that the possibility of a strongly cor- 
related Kondo ground state possessing an SU(4) sym- 
metry has also been discussed very recently in the dif- 
ferent context of two small dots coupled with a strong 
capacitive inter-dot coupling^ The possibility of orbital 
and spin Kondo effects in such a geometry was previ- 
ously anticipated by Schon et alwL inspired by prelim- 
inary experiments of Ref. (2^]. It is worth noting that 
these types of problems have also potential connections 
with the twofold orbitally degenerate Anderson impurity 
model 30 ' 31 and more precisely with the physics of cer- 
tain heavy fermions like UBei3 where the U ion is mod- 
eled by a non-magnetic quadrupolar doublet^ and then 
quadrupolar (orbital) and spin Kondo effects can in prin- 
ciple interfered 

Our paper is structured as follows: In Section II, we 
resort to a Schrieffer- Wolff transformation and derive the 
effective model including the different useful parameters. 
In Section III, assuming that we are far from the degener- 
acy points of the grain we use a pedestrian perturbation 
theory; This reveals the importance of spin flips even in 
this limit. In Section IV, we carefully investigate both 
theoretically and numerically the interplay between or- 
bital and spin Kondo effects at the degeneracy points. 
In Section V, we discuss in details the effects of possi- 
ble symmetry breaking perturbations and the crossovers 
generated by such perturbations. Finally, section VI is 
devoted to the discussion of our results and especially we 
summarize our main experimental predictions for such a 
setup. 
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II. MODEL AND SCHRIEFFER- WOLFF 
TRANSFORMATION 

In the sequel, we analyze in details the behavior of 
charge fluctuations in the grain. In order to model the 
setup depicted in Figure 1, we consider the Anderson- like 
Hamiltonian: 



H 



Q 2 



k p 
+ ta\a„ + Un-\ny (1) 

kcr pa 

where etfc CT , a CT , a pa are the annihilation operators for 
electrons of spin a in the lead, the small dot, and the 
grain, respectively, and t is the tunneling matrix element 
which we assume to be k independent for simplification. 
Let us first consider that tunnel junctions are symmet- 
ric. We also assume that junctions are narrow enough 
and contain one transverse channel only. Extensions of 
the model to asymmetric or larger junctions will be an- 
alyzed later in section V. We also assume that the en- 
ergy spectrum in the grain is continuous, which implies 
that the grain is large enough such that its level spac- 
ing A s is very small compared to its charging energy 
E c = e 2 /(2C): A g /E c -► (in Ref. 0A 9 ~ 70mK was 
not sufficiently small compared to the Kondo temper- 
ature scale which hindered the logarithmic capacitance 
peakS to completely develop) . Q denotes the charge oper- 
ator of the grain, C is the capacitance between the grain 
and the gate electrode, and ip is related to the back-gate 
voltage V g through ip = — V g . e < and U are respec- 
tively the energy level and charging energy of the small 



dot, and 



= «t, 



The inter-dot capacitive coupling is 



assumed to be weak and therefore neglected. 

We mainly focus on the particularly interesting situa- 
tion where the small dot is in the Kondo regime, which 
requires the last level to be singly occupied and the con- 
dition 



t < 



(2) 



to be satisfied (e < 0). The resonant level limit where e 
lies near the Fermi level will be addressed at some points 
in Section V. In the local moment regime, we can in- 
tegrate out charge fluctuations in the small dot using 
a generalized Schrieffer- Wolff transformation.iSiS More 
precisely, the system is described by the Hamiltonian: 

Q 2 

H = e*4 * + e P a P a P + ^ + ¥Q (3) 



■ ^ + V ) a ™a«- 



To simplify the notations, the spin indices have been 
omitted and hereafter, m, n take values in the two sets 



"lead" (k) or "grain" (p), the spin S is the spin of the 
small dot, a are Pauli matrices acting on the spin space 
of the electrons. Let us now discuss the parameters J 
and V in more details. 

In the vicinity of one degeneracy point obtained for 
ip = — e/2C, where the grain charging states with Q = 
and Q — e are degenerate, we find explicitly: 



J = 2t 2 



1 



1 



U 



A small direct hopping term 



V=1 2 



1 



U 



(4) 



(•5) 



is also present and should not be neglected. In particular, 
this embodies the so-called "charge flips" from the reser- 
voir to the grain and vice- versa in the Matveev's original 
problem. Notice that the ratio V/J can take values be- 
tween -1/4 (when U = —e) and 1/4 (when U — > oo). 
V = corresponds to the particle-hole symmetric case 
where 2e + U = 0. For ip = — e/2C, the energy to add a 
hole or an electron onto the metallic grain vanishes and 
therefore the Schrieffer- Wolff parameters V and J are 
completely identical to those of a small dot connected to 
two metallic reservoirs. 34 Furthermore, remember that in 
the present model the ultraviolet cutoff at which the ef- 
fective model becomes valid can be roughly identified to 
D ~ minli^c, A^} where is the level spacing of the 
small dot (with today's technology it is possible to reach 5 
A d - 2-3K and for the gram 2 Q E c - 2.3K). 

On the other hand, far from the degeneracy point tp = 
—e/2C - which means on a charge plateau - the energy 
to add a hole on the grain is U-i — E c (l + 2N) where 
N = CVg/e 1/2. Similarly, it costs Ui = E c (l - 2N) 
to add an extra electron onto the grain. The lead-dot and 
grain-dot Kondo couplings, Jo and J\ respectively, then 
become asymmetric even for symmetric junctions: 



Jo = 2t z 



Ji = 2t z 



1 



1 



-e U + e 
1 



J 



(6) 



Ux-t U + e + U- 



In the second equation, the virtual intermediate state 
where an electron first hops from the grain onto the small 
dot induces an excess of energy U—x in the second term. 
The first term contains the energy of the intermediate 
state of the process where the temporal order of the hop- 
ping events is reversed. The off-diagonal terms where an 
electron from the reservoir [grain] flips the impurity spin 
and then jumps onto the grain [reservoir] reads 



Joi — 2t 



Jio — 2< 



1 



1 



U t -t U + e 
1 1 



(7) 



-e U + e + U- 
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Note that in general particle-hole symmetry is absent in 
the large dot, so in principle Joi ^ Jio- But, in our set- 
ting, E c = e 2 /2C <C |e|, U + e, so in the following we will 
neglect the asymmetry between Joi and Jio far from the 
degeneracy points (Jqi = Jio) (this has no drastic conse- 
quence on the results) . In the finite temperature range 
T < U\, U-i, these off-diagonal processes are suppressed 
exponentially as Jio = Jxo(T) s» Jioe~ Ul / ikT , whereas 
the diagonal spin processes can be strongly renormal- 
ized at low temperatures. In other words, in the Renor- 
malization Group language, if we start at high temper- 
ature with a set of Kondo couplings Jo, Ji, Joi, JiO) the 
growing of Joi, Jio is cut-off when T is decreased below 
max([/i, U—i) whereas the growing of Jo, Ji is not. This 
offers a room to reach a 2-channel Kondo effect in the spin 
sector (for asymmetric tunneling junctions) provided the 
condition J = Ji can be reached with a fine-tuning of 
the gate voltages^ We can do the same approximation 
for the V term and define Vio, Vbi accordingly (with 
Vio — Voi) and also Vio. 



III. PEDESTRIAN PERTURBATION THEORY 
ON A PLATEAU 

We want first to compute the corrections to the aver- 
age charge on the grain on a charge plateau due to the 
Kondo and V couplings bearing in mind that when the 
tunneling amplitude t — > 0, the average grain charge (Q) 
exhibits perfect Coulomb staircase behavior as a func- 
tion of V g . We confine ourselves to values of tp in the 
range — e/(2C) < tp < e/(2C), which corresponds to the 
unperturbed (charge) value Q = 0. A first natural ap- 
proach is to assume that the Kondo and charge-flip cou- 
plings are very small compared to the charging energy 
E c = e 2 /(2C) of the grain and to calculate the correc- 
tions to Q = in perturbation theory. Despite this per- 
turbative calculation will appear of limited use, it is very 
instructive to perform it in order to point the different 
sources of divergences that appear when approaching the 
degeneracy points, the main issue treated in this paper. 
At second order, we find 



(Q) a =e -J 



'10 



2V 



10 



Note that at finite low temperature T < U\,U-\, 
we should use the renormalized off-diagonal couplings 
■Aoi Vio which are small (in other words the flow of 
the off-diagonal Kondo couplings has been cut-off for 
T < Ui,U-i). This better reproduces the (exact) nu- 
merical calculations of Ref . For more details, we refer 
the reader to Appendix A. The density of states in the 
lead and in the grain have been assumed to be equal 3 
and taken to be 1 for simplicity. This result tends to triv- 
ially generalize that of a grain directly coupled to a lead. 9 
However, there are two reasons that may suggest this per- 
turbative approach is divergent. Higher-order terms - al- 
ready at cubic order- involve logarithmic divergences as- 



sociated to the renormalizations of the Kondo couplings 
(see Appendix A), but also other logarithms indicating 
the vicinity of the degeneracy point tp = — e/2C in the 
charge sector. For example, a correction at cubic order 
to the result in Eq. © is given by 



(Q) 3 oc J Jio 2 In 



D 



In 



e/2C-<p 
e/2C + tp 



(9) 



We also have a similar correction in JiJio 2 . It would 
be potentially interesting to observe the logarithmic 
temperature-dependence of (Q) on a given plateau due to 
Kondo spin-flip events. Note also that the perturbation 
theory in the Vio term has been previously extended to 
the fourth order. 10 The perturbative result is valid only 
far from the degeneracy points provided the renormaliza- 
tion, e.g., of the spin Kondo coupling Jo is also cut-off 
either by the temperature T or by a magnetic field B [in 
general, for symmetric junctions one already gets Jo > Ji 
at the bare level; See Eq. (6)]. This considerably restricts 
the range of application of this perturbative calculation 
compared for example to the simpler setup involving a 
grain coupled to a reservoir and even on a charge plateau 
the temperature much be larger than the emerging spin 
Kondo energy scale between the lead and the small dot. 
Finally, note that in our perturbative treatment at fi- 
nite temperature T < U±, £7_i we have made the (stan- 
dard) approximation: We have only virtually introduced 
the temperature through the renormalization of the cou- 
plings Jio and Vio- 

The other regimes which requires non-perturbative ap- 
proaches will be studied in Sections IV and V. 



IV. ORBITAL AND SPIN KONDO EFFECTS 
CLOSE TO THE DEGENERACY POINTS 

In this section, we will be primarily interested in the 
situation close to the degeneracy point tp — — e/2C where 
none of the perturbative arguments above can be applied. 
We want to show that the Hamiltonian given by Eq. © 
can be mapped onto some generalized Kondo Hamilto- 
nian following Ref. 



A. Mapping to a generalized Kondo model 

Close to the degeneracy point tp = —e/2C and for 
ksT -C E c , only the states with Q — and Q — e are 
accessible and higher energy states can be removed from 
our theory introducing the projectors Po and P\ (which 
project on the states with Q = and Q — e in the grain 
respectively). The truncated Hamiltonian © then reads: 



H = 



£ka kT a kT 

fc,T=0,l 



P + Pi)+ ehPi 



(10) 



k,k' 



2 a 



S + V 



(4 



ak'oPo + al, a k iPi 



5 



E 

T=0,1 



J 



a ■ S + V I a\ T a k 



where now the index r = indicates the reservoir and 
r = 1 indicates the grain. We have also introduced the 
small parameter 

which measures deviations from the degeneracy point. 
Considering t as an abstract orbital index, the Hamil- 
tonian can be rewritten in a more convenient way by 
introducing another set of Pauli matrices for the orbital 
sector&ii 



H 



£ k a kT a kj 



E 

k,k' 



+ 



E 

T,T 

E 



j. 



S + V) ( t x T x + t v T v 



(12) 



')< 



TjT ,a' kT a k >i 



J 



3 ■ S + V ) a kr a k 



In this equation, the operators (S, a) act on spin and the 
(T, r) act on the (charge) orbital degrees of freedom. 

The key role of this mapping stems from the fact that 
(Q) can be identified as (an orbital pseudo-spin) 



(Q) 



(T z ) 



(13) 



Then, we can introduce the extra (charge) state \Q) as an 
auxiliary label to the state |$) of the grain. In addition 
to introducing the label \Q) we make the replacement 



a fci a fc'o-Po 

a k O a k'lP\ 



a kl a k >oT^ 
a k \T~ 



(14) 



b k'0 



Notice that T + and T~ are pseudo-spin ladder operators 
acting only on the charge part \Q). More precisely, we 
have the correct identifications 



T~\Q 
T+\Q 



1) 
0) 



T-\T Z 
T+\T Z 



+1/2) 
-1/2) 



\Q 
\Q 



0} (15) 
1), 



meaning that the charge on the single- electron box is 
adjusted whenever a tunneling process takes place. Fur- 
thermore, since T+\Q = 1) = and T~\Q = 0) =0 
these operators ensure in the same way as the projection 
operators Po and P\ that only transitions between states 
with Q = and Q — 1 take place. This leads us to 
identify Pi + Pq with the identity operator on the space 
spanned by |0) and |1) and P\ — Pq with 2T Z . We now 
introduce an additional pseudo-spin operator via: 



t 1 t - 

a kl a k >Q = - a' kT T a k , T , 

t 1 t + 

a kO a k'l = 2 a fcr T a k'T>, 



(16) 



where the matrices = t x ± ir y are standard com- 
binations of Pauli matrices. Finally, the Coulomb 
term h mimics a magnetic field acting on the orbital 
space. Therefore, the (quantum) grain capacitance C q — 
~d(Q)/dh is equivalent to the local isospin susceptibil- 
ity XT = —d(T z )/dh up to a factor e. For simplicity, 
we will substract the classical contribution C which is 
V^-independent . 

But obviously, to compute the latter, we have to deter- 
mine the nature of the Kondo ground state exactly. 

Typically, when only "charge flips" are involved 
through the V term, the model can be mapped onto a 
two-channel Kondo model (the two channels correspond 
to the two spin states of an electron), and the capac- 
itance always exhibits a logarithmic divergence at zero 
temperature^ Here, we have a combination of spin and 
charge flips. Can we then expect two distinct energy 
scales for the spin and orbital sectors? To answer this 
question, we perform a perturbative scaling analysis fol- 
lowing that of a related model in Ref. 36]. We first 
rewrite the interacting part of the Hamiltonian in real 
space as: 

+ Y±T z (^T z ^) + y±[T+(^T-^)+h.c.] (17) 
+ Q Z T Z S ■ {^T z a^) + Q^S ■ [T+i^T-ffifi) + h.c] , 

where V> T(T = J2k a krcr- 

A host of spin-exchange C3> isospin-exchange interac- 
tions are clearly generated; J refers to pure spin-flip pro- 
cesses involving the S=l/2 spin of the small dot, V± to 
pure charge flips which modify the grain charge, and Q± 
describes exotic spin-flip assisted tunneling. 

This Hamiltonian exhibits a structure which is very 
similar to the one introduced in Ref. |26| in order to 
study a symmetrical double (small) quantum dot struc- 
ture with strong capacitive coupling. 26 However, since 
the physical situation that led us to this Hamiltonian 
here is very different from that of Ref. 26], our bare 
values for the coupling parameters are also very different 
(forJ<l): 

V± = V, V z = , Q z = , Q± = J/4. (18) 

We have ignored the potential scattering Vip^ip which 
does not renormalize. It is also relevant to note that 
this model belongs to the general class of problems of 
two coupled Kondo impurities. However the coupling 
between impurities, namely Q±, is far different from the 
more usual RKKY interaction^ 

Again, bear in mind that here the operators Pi.o = 
(1 ± 2T z )/2 and p ,i = (1 ± t z )/2 project out the grain 
state with Q = e and Q = 0, and the reservoir/grain 
electron channels, respectively. The spin S corresponds 
to the spin of the small dot in the Kondo regime and the 
index a is the spin state of an electron in the reservoir or 
in the grain. 
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Note that in the situation of Ref. 26], the operators 
P± = (l± 2T Z ) /2 and p± = (1 ± t z )/2 rather project out 
the small double dot states (n+,n_) = (1,0) and (0,1), 
and the right/left (+/— ) lead channels, respectively. Ad- 
ditionally, the spin 5* is the spin (excess) impurity either 
on the left or the right dot and the index a denotes the 
spin state of electrons in the reservoirs. The correspond- 
ing bare values in that case would be rather of the form: 



V± = Qj 



V z 



= J. 



(19) 



B. Perturbative Renormalization Group analysis 

The low-energy Hamiltonian can be treated using per- 
turbative renormalization group (RG). It is relevant to 
observe that no new interaction terms are generated to 
second order as the bandwidth is reduced. By integrating 
out conduction electrons with energy larger than a scale 
E -C D (~ min{.E c , Ad} being either the level spacing 
of the small dot or the charging energy of the grain, i.e., 
the ultraviolet cutoff) , we obtain at second order the fol- 
lowing RG equations for the five dimensionless coupling 
constants: 



dJ 

dl 
dV x 
~dT 
dQ z 

dl 
dQ±_ 

dl 



J 2 + Q 2 + 2Q ± 2 



V ± 2 + 3Q± 2 
V X V Z + 3Q X Q Z 
2JQ Z + 2V ± Q ± 
2JQ± + V Z Q ± + V ± Q Z , 



(20) 



with I = ln[D/E] being the scaling variable, E is the 
running bandwidth. This RG analysis is applicable only 
very close to the degeneracy point ip — — e/(2C) where 
the effective Coulomb energy in the grain or h vanishes 
and obviously only when all coupling constants stay <C 1. 
Higher orders in the RG have been neglected. 

Although the equations l|20|l have no simple analytic 
solution, one can try to read off the essential physics from 
numerical integration and the initial conditions (|18|) . 

Let us first discuss the most obvious case of a particle- 
asymmetric level, with V± > meaning (large) U ^> — 2e. 
In this case, the numerical integration of the RG flow in- 
dicates that even though we start with completely asym- 
metric bare values of the coupling constants, all cou- 
plings diverge at the same energy scale due to the pres- 
ence of the spin-flip assisted tunneling terms Q± and Q z . 
This energy scale that we can identify with a generalized 
Kondo temperature is difficult to calculate analytically. 
However, we can approximate it by the one of the com- 
pletely symmetrical model 



Furthermore, we have checked numerically that all cou- 
pling ratios converge to one in the low energy limit 
provided the RG equations can be extrapolated in this 
regime. These results have been summarized in Figure 2. 
As confirmed below with an exact Numerical RG treat- 
ment, the entanglement of spin and orbital degrees of 
freedom in this geometry will lead to an higher symme- 
try than SU(2) ® SU(2), namely SU(4), and then to the 
formation of a Fermi-liquid correlated ground-state with, 
e.g., the complete screening of the orbital spin T. [SU(4) 
is the minimal group allowing spin-orbital entanglement 
and which respects rotational invariance both in spin and 
orbital spaces.] Recall that the presence of the spin-flip 
assisted tunneling terms then definitely hinders the pos- 
sibility of a non- Fermi liquid ground state induced by the 
over-screening of the the pseudo-impurity T. 
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2000 2500 3000 3500 

FIG. 2: Evolution of the four coupling ratios as a function 
of the scaling variable I — \n(D/E). The initial conditions 
have been chosen as: J(0) = u, Vj_(0) = O.lOu, Q±(0) = u/4 
with u = 0.00018 and Vi(0) = Q z (0) = 0. The full line is 
Q±/ J, the dotted line Vj_/V z , and the dashed lines to Q±/V± 
and Q_l/Q z (which diverges for I — » 0). All the couplings 
are strongly renormalized for l c m 3914 and all their ratios 
converge to 1. Extrapolating the flow to I ^> l c would give a 
straight horizontal line where the coupling ratios remain 1. 

Let us now analyze the particle-hole symmetric case, 
i.e., V± — 0. At second order, the RG flow would tend to 
suggest that two parameters, namely V± and Q z remain 
zero whatever the energy scale. Typically, the Kondo 
coupling J is the largest throughout the RG flow and 
seems to be the first one to diverge. On the other hand, 
the ratios V z /J and Q±/J cannot be neglected which 
tends to exclude an SU(2) x 5/7(2) symmetry where the 
spin and the orbital degrees of freedom would be inde- 
pendently screened (Figure 3). Instead, spin-orbital mix- 
ing (entanglement) seems to be prominent at low-energy. 
Even though the perturbative RG is certainly not suffi- 
cient to draw more definitive conclusions, it is also in- 
structive to observe that for V± < 0, the ratios Q±/J 
and V z /J still converge to one. Since the system defi- 
nitely has to restore the rotational invariance both in spin 
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FIG. 3: Here, we have chosen: J(0) = u, Q±(0) = u/4 with 
u = 0.00018 and Vj_(0) = 14(0) = Q z (0) = 0. The coupling 
J(l) is the largest throughout the RG flow, but the ratios 
Q±/J and V z /J cannot be neglected. Furthermore, at sec- 
ond order, the couplings V± and Q z would remain zero. How- 
ever, the NRG concludes that even in this limit the system 
converges to an SU(4) Fermi- liquid fixed point with identi- 
cal coupling constants which emphasizes the importance of 
higher-order terms and that spin-orbital mixing is well promi- 
nent and that the rotational invariance is restored both in spin 
and orbital spaces. 



and orbital spaces, this tends to emphasize that higher- 
order terms play a crucial role in the crossover regime 
eventually by restoring an SU(4) Fermi-liquid even for 
those cases. Moreover, the RG analysis suggests that 
the temperature scale at which the Fermi-liquid behav- 
ior emerge would be much smaller for vanishing and neg- 
ative V± because the system needs a much longer time 
to restore the rotational symmetry both in spin and or- 
bital spaces. To enumerate higher order terms would 
be a very tedious task, therefore this assertion will be 
rather checked by NRG, a completely non perturbative 
method. To summarize this part, we emphasize that for 
V±_ < 0, the above perturbative analysis does not allow us 
to determine the precise nature of the low-temperature 
fixed point, whether the orbital (isospin) moment is ex- 
actly screened or over-screened. We will prove in sec- 
tion ^^E] using NRG, that a SU(4) strongly-correlated 
ground state emerges for any physical value of Vj_, i.e., 
-J/4 <V±< J/4. 



Entanglement of spin and charge degrees of 
freedom 



This RG analysis suggests - at least for not too small 
positive V± - that our model becomes equivalent at low 
energy to an SU(4) symmetrical exchange model: 



a/3 



T f3 



(22) 



4L M 
4Z T 



A \ " 



A 



Since all the coupling ratios converge to one, we have 
rewritten the Kondo Hamiltonian l|17fl with the unique 
coupling constant J. We have introduced the "hyper- 
spin" 

M A e {2S a ,2T a ,4S a T> 3 } , (23) 

for a,/3 = x,y,z. The operators M A can be regarded 
as the 15 generators of the SU(4) group. Moreover, this 
conclusion will be strongly reinforced by the NRG analy- 
sis proposed below (whose range of validity is broader 
than Eqs. pOjl l which indeed concludes that the ef- 
fective Hamiltonian (|22|l is appropriate for all values of 
-J/4 < V± < J/4. Note that apparently has 
SU(2) x SU(2) symmetry, representing rotational in- 
variance in both spin and orbital (pseudo-spin) spaces, 
and also interchange symmetry between spin and pseudo- 
spin. But, the full symmetry is actually the higher sym- 
metry group SU(4), which clearly unifies (entangles) the 
spin of the small dot and the charge degrees of freedom 
of the metallic grain. Notice that the irreducible repre- 
sentation of SU(4) written in Eq. IL'l'l) has been used 
previously for spin systems with orbital degeneracy^*^ 
The electron operator ip now transforms under the fun- 
damental representation of the SU(4) group, with gener- 
ators t A (A = 1, 15), and the index \x labels the four 
combinations of possible spin (T,|) and orbital indices 
(0, 1), which means (0, T), (0, |), (1, T) and (1, |). 

The emergence of such a strongly-correlated SU(4) 
ground state, characterized by the quenched hyper-spin 
operator 



* + i f + i 



(24) 



clearly reflects the strong entanglement between the 
charge degrees of freedom of the grain and the spin de- 
grees of freedom of the small dot at low energy induced 
by the prominence of spin-flip assisted tunneling. There 
is the formation of an SU(4) Kondo singlet which is a 
singlet of the spin operator, the orbital operator, and the 
orbital-spin mixing operator U a ^ — S a T 13 . Again, let us 
argue that this enlarged symmetry arises whatever the 
parameter V± simply because the spin-flip assisted tun- 
neling term Q± always flows off to strong couplings at 
the same time than the more usual Kondo term J; The 
system then must inevitably converge to a fixed point 
with orbital-spin mixing. To respect rotational invari- 
ance in both spin and orbital spaces the only possibility 
is indeed an SU(4)-symmetric Kondo model (as agreed 
with NRG). 



D. Capacitance: Destruction of Matveev's 
logarithmic singularity 

The (one-channel) SU(N) Kondo model has been ex- 
tensively studied in the literature (see, e.g., Ref. jUj). 
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In particular, the strong coupling regime corresponds to 
a dominant Fermi liquid fixed point induced by the com- 
plete screening of the hyper-spin M a , implying that all 
the generators of SU(4) yield a local susceptibility with a 

behavior in 41 . ~ l/T^F^. T z being one of these gener- 
ators, we deduce that \t — —d(T z }/dh and then the 
(quantum) capacitance of the grain C q = —d(Q)/dh 

roughly evolves as 1/T^ U ^ at low temperatures^ We 
have substracted the classical capacitance C. Conse- 
quently, for h <C e/C, we obtain a linear dependence 
of the average grain charge as a function of V g — —tp: 

h e 



(Q)-- 2 = 



SU(4) 



K 



T 



SU(i) 
K 



e 
2C 



(25) 



The hallmark of the formation of the SU(4) Fermi liq- 
uid in our setup is now clear. The (grain) capacitance 
peaks are completely smeared out by the mixing of spin 
and charge flips and Matveev's logarithmic singularity 9 
has been completely destroyed. Additionally, the strong 
renormalization of the V (and J) term - and the stabil- 
ity of the strong-coupling Kondo fixed point - clearly re- 
flects that the effective transmission coefficient between 
the lead and the grain becomes maximal close to the 
Fermi level. (The maximum of the tunneling appears 
not exactly at the Fermi level as one could guess from 
the value of the phase shifts S = 

This example could also be interpreted as an interest- 
ing proof that one can already wash out the Coulomb 
staircase when the 'effective' transmission coefficient be- 
tween the grain and the lead is roughly one only close to 
the Fermi energy (and not for all energies^). Conceptu- 
ally, this is not accessible with a small dot in the resonant 
level limit ^2*24 We stress that this stands for a remark- 
able signature of the formation of a Fermi-liquid ground 
state when tunneling through a single- electron box. 



E. Confirmation with Numerical Renormalization 
Group analysis 

In order to confirm the results obtained by pertur- 
bative RG and extend our investigation to the strong 
coupling regime, we have performed a collaborative 
NRG 33,42 analysis of the model described by Eq. lfT7|) 
similar as in Ref. [2(j. Note in passing that the model 
of Eq. (|17|l with asymmetric bare values is not strictly 
speaking integrable. Therefore, we resort to the NRG 
method which in general can be successfully applied to 
(various) two-impurity Kondo models^ 3 . At the heart of 
the NRG approach is a logarithmic energy discretization 
of the conduction band around the Fermi points. In this 
method - after the logarithmic discretization of the con- 
duction band and a Lanczos transformation - one de- 
fines a sequence of discretized Hamiltonians, Hat, with 
the relation* 4 ^ 



H 



N+l 



N 



JV 



(/iV,TO-/w+l 



h.c 



(26) 



where / , TCT = 4> T a/V2 and H = 2A 1 / 2 /( 1 + A ) h k with 
A ~ 3 as discretization parameter, and £jv ~ 1- For 
the definition of /at see Ref. ji^. The original Hamilto- 
nian is connected to the -Hat's as H — Yluin^oo ^nHn 
with lu n = A- (N+1 ^ 2 (1 + A)/2. Using the logarithmic 
separation of the energy scales we are allowed to diago- 
nalize -Hat's iteratively and calculate physical quantities 
directly at the energy scale u> ~ uin ■ We have calculated 
the dynamical spin and orbital spin (ac) susceptibilities 



9fmxoH = 9fm^<[O(t) J O(0)]) 



(27) 



where O = T z , S z and J- denotes the Fourier transform. 
According to the discussion above, the couplings were 
chosen as J = 4Q^, Q z = V z = 0. 

The obtained orbital spin susceptibility for different 
values of V± is shown in Figure^] Regardless of the value 
of V±, the T z susceptibility exhibits a typical Fermi- 
liquid like peak at an energy scale which can be identified 

as Tx U{4) . Above this energy scale it behaves as \ ~ to 
indicating that the correlation function in Eq. I|27|) is con- 
stant for very short times while for u> < T^- 11 ^, \ ~ 10 
as a signature of the ~ 1/t 2 asymptotic of the afore men- 
tioned correlation function for a Fermi-liquid model. In- 
deed, at T=0, this ensures a grain's capacitance, 



C q = 



1/T SU(4) 



dt ([T z (t),T z (0)}) = 



1 



T, 



5(7(4) 



(28) 



A" 



Furthermore, as one can see in Figures 01 (for A 2 — > 0) 
the Kondo screening simultaneously takes place in the 
spin and orbital sectors, indicating the S'H(4)-symmetric 
nature of the effective low energy Hamiltonian. 




FIG. 4: The orbital spin T z susceptibility for different values 
of V_l . In all the cases the susceptibility shows a typical SU(4) 

Fermi-liquid state at lu = T^ U ^ (V±). Inset : As a compari- 
son we plot the same quantity for the 2-channel Kondo model. 



Furthermore, we can clearly observe that T 1 



SU{4) 



markedly 



decreases for lower values of U meaning V±/J < (i.e., by 
making the small dot larger and larger 4-) 
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To give a rigorous proof of the SU (4) Fermi liquid ground 
state one has to analyze the finite size spectrum obtained 
by NRG. It turns out that (as in Ref. |2f|) the spectrum 
can be understood as a sum of four independent chiral 
fermion spectra with phase shift 7r/4 in accordance with 
the prediction of the SU(A) Fermi liquid theory. This 
result proves that the low-energy behavior is described 
by the Fermi liquid theory even at V± = 0, but as con- 
jectured above the temperature scale at which the Fermi 
liquid emerges decreases as we change the coupling V± 
from 0.4J to -0.4J. 

For comparison, in the inset of Figurc^Jwe plot the dy- 
namical susceptibility for the two channel Kondo model: 
In that case, 9?mx(w) ~ const, which in contrast traduces 
that the capacitance C q would exhibit a logarithmic di- 
vergence at zero temperature^ 

Additionally, the SU(4) Kondo temperature scale is 
considerably reduced for negative values of Vj_, i.e., by 
decreasing the on-site interaction U on the small dot 
(U <C — 2e). This makes sense since by substantially 
decreasing the Coulomb energy of the small dot - i.e., 
by progressively increasing the size of the small dot - 
one expects the breakdown of the SU(4) fixed point and 
a situation similar to that of a reservoir and two large 
dots^i (According to Eq. J2J), spin Kondo physics should 
definitely vanish for [/ <SC — e). 



ON THE STABILITY OF THE SU(4) FIXED 
POINT AND CROSSOVERS 



In contrast to the two-channel Kondo fixed point, 
which is known to be extremely fragile with respect to 
perturbations (e.g., channel asymmetry, magnetic field), 
the SU(4) fixed point is robust at least for weak pertur- 
bations. 

In order to demonstrate the robustness of the SU(4) 
Fermi liquid fixed point we have checked the role, e.g., of 
a magnetic field in real and orbital spin sectors. It turns 
out that both terms are marginal operators in RG sense. 
On the other hand, when the magnetic [orbital] field is 
much larger than the scale of the Kondo temperature, 
the processes which involve spin [orbital spin] flips are 
suppressed and low energy physics is described by a one- 
channel orbital spin [spin] Kondo effect, with a smaller 
Kondo temperature than that of the SU(4) case. Let 
us now thoroughly analyze the different fixed points, the 
effects of an asymmetry between the tunnel junctions and 
of rather large junctions with more conducting channels. 



A. Magnetic field 

First of all, we have checked with NRG that the SU(4) 
Fermi liquid fixed point resists for quite weak external 
magnetic field. But, applying a strong magnetic field 
B ^> Tk unavoidably destroys the SU(4) symmetry. But 
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FIG. 5: The orbital spin T z susceptibility for different val- 
ues of the external magnetic field B. The low-energy physics 
consists of a Fermi liquid regardless of B, but the symmetry 
is reduced for large magnetic fields to SU(2) (for the orbital 
space) and the Kondo energy scale as well. 



at zero temperature, we expect the behavior of charge 
fluctuations close to the degeneracy points to remain 
qualitatively similar. Indeed, in a large magnetic field 
spin flips are suppressed at low temperatures, i.e., Q± = 
Q z = J = 0, and the orbital degrees of freedom, through 
V± and V z , develop a standard one-channel Kondo model 
(the electrons have only spin- up or spin down), which also 
results in a Fermi-liquid ground state with a linear de- 
pendence of the average grain charge as in Eq. Q25[l. Yet, 
the emerging Kondo temperature will be much smaller, 



T K [B 



(29) 



with for instance V ~ t 2 /(—2e) for U — > +oo, and might 
not be detectable experimentally. A substantial decrease 
of the Kondo temperature when applying an external 
magnetic field B has also been certified using NRG even 
for extremely large values of V (Figure 5). 



B. Away from the degeneracy points: small dot as 
a resonant level 



A weak orbital magnetic field (orbital splitting) A z oc 
h does not modify the SU(4) Fermi liquid state. 

Moreover, the application of a strong A z always leads 
to a single-channel Kondo effect in the spin sector. A 
naive consideration - focusing on the RG flow above 
- would suggest the possibility of a two-channel (spin) 
Kondo effect: The simultaneous screening of the excess 
spin of the small dot by the lead and the grain electrons, 
independently. However, going back to the Schrieffer- 
Wolff transformation for the situation away from the de- 
generacy points, the charging energy of the metallic grain 
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FIG. 6: The real spin S z susceptibility for different values 
of the orbital splitting, A z . For A z > T^ U{4) the processes 
which involve orbital spin flip are suppressed resulting in a 
purely one-channel spin Kondo effect with a smaller Kondo 
temperature of the order of that for a small dot embedded 
between two leads 2V[A„]. Recall that the energy scale at 
which the SU(4) correlated state arises can be much larger 
than Tff[A z ] which should certainly ensure the observation of 
our theoretical results. It is worthwhile to note the parallel 
between Figures 5 and 6 by interchanging T z <-> S z and A z <-> 
B (however, T K [A Z ] > T K [B]). 

definitely ensures J\ ^ Jo (provided we start with al- 
most symmetric junctions), a condition that destroys the 
stability of the two-channel spin Kondo fixed point. The 
spin Kondo coupling J will be the first one to flow off to 
strong couplings (as anticipated in Section III) . The NRG 
calculation clearly confirms this expectation: the A z 
term not only suppresses the orbital spin-flip terms but 
also generates an asymmetry between the grain-dot and 
lead-dot spin couplings which destroys the two-channel 
Kondo behavior. The possible two-channel (spin) Kondo 
regime proposed by Oreg and Goldhaber-Gordon2& can 
not be reached with this model, at least, for symmetric 
junctions. Asymmetric junctions and a fine tuning of the 
grain gate voltage far from the degeneracy points would 
be necessary to reach the condition J = J\. On the 
other hand we will see that, for quite asymmetric barri- 
ers, a two-channel Kondo behavior rather for the orbital 
degrees of freedom can appear near the degeneracy points 
but at extremely small (and a priori unreachable) tem- 
peratures. 

For A z ^> T^- U< " 4 \ the Kondo temperature scale here 
resembles that for a small dot connected to two leads 
(Jo = J)2A and, in principle, is still experimentally ac- 
cessible: 

T K [A Z ] ~ D e~ ljJ < T| a(4) . (30) 

Henceforth, this will cutoff the logarithmic divergence in 
the charge fluctuations away from the degeneracy point 



FIG. 7: The orbital spin T z susceptibility for different val- 
ues of the orbital splitting, A z . For A z > T^ U ^ the pro- 
cesses which involve orbital spin flip are clearly suppressed 
at a scale of A z producing instead a Schottky anomaly. The 
orbital pseudospin model then becomes inappropriate to de- 
scribe the charge fluctuations of the grain at low energy. We 
rather apply another resonant level mapping and a perturba- 
tion theory similar to that of Ref. [24|. 

ip = —e/2C [see Eq. (9)]. In order to describe the physics 
at strong orbital magnetic field, i.e., away from the de- 
generacy points and at lower temperature and more pre- 
cisely the average grain charge (Q) , we seek to go beyond 
the effective model in Eq. I|17fl . Indeed, at energy smaller 
than Tk[A z ], the physics can be qualitatively identified 
with that of Ref. [24j : The Kondo screening of the excess 
spin of the small dot by the lead produces an Abrikosov- 
Suhl resonance at the Fermi level, and the small dot plus 
the lead can be replaced by a resonant level with the en- 
ergy e — ► and the resonance width ~ Tk[A z ). Now, one 
can still allow for a (weak) residual tunneling matrix ele- 
ment t between the grain and the effective resonant level 
(which may be of the same order as the bare tunneling 
matrix element t between the small dot and the grain but 
its value is difficult to determine accurately). For an illus- 
tration, see Figure 8. Reformulating results of Ref. [24j 
for our case and including that Tk[A z ] <C U\,U-i for 
N = CVg/e < 1/2 (ip < -e/2C), at zero temperature 
we find 

(31) 

r an 

6 E c ir (1-2A)(1 + 2AT)' 
with the effective tunneling energy scale 

r = 7r5^t 2 <y(e p )«:l7i,E7_i. (32) 
p 

Since Ui and U-i are of the order of E c for N <C 1/2, we 
observe that the charge smearing far from the degeneracy 
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FIG. 8: Illustrative view of the effective low-energy model for 
almost symmetric barriers away from the degeneracy points: 
According to Eq. JSJ the charging energy on the grain in- 
evitably ensures that the spin Kondo coupling Jo between the 
bulk lead and the small dot will be the first one to flow off to 
strong couplings at the energy scale Tfr[A z ]. The grain be- 
comes virtually weakly-coupled to an effective resonant level 
with a reduced bandwidth ~ Tk[A z ] <C D. 

points is small at low temperatures. Additionally, recall 
that for JV< 1/2 and zero temperature, at second order 
in t the average grain charge also exhibits a (small) linear 
behavior as a function of N or V g which is slightly distinct 
from the original Matveev's situation (Figure 9)AiS 



SU(4) Fermi liquid 
linear slope also for N=l/2 






2 







1 


^ y 


-0.3 -0 J .2^J£J^ > 
. — ^ — 
/ 


1 


0.1 0.2 0.3 


/ "° 


2 





effective resonant level model 
small linear slope for N=0 



FIG. 9: Profile of the average charge (Q) on the grain 
versus N — CV g /e for (almost) symmetric junctions and 
T < Tk[A z ]. Again, the SU(4) Kondo entanglement between 
spin and orbital degrees of freedom, e.g., at the degeneracy 
point N = 1/2 produces a Fermi-liquid state and the Coulomb 
staircase exhibits a conspicuous smearing. Away from the 
degeneracy points, the physics becomes similar to that of a 
resonant level weakly coupled to a grain which also ensures a 
linear (but small) behavior for [Q)[N] when N -> 0; The full 
line curve corresponds to T/E c = 0.15 and the dashed line 
curve to T/E c = 0.1. 



C. Case of asymmetric junctions 

Another interesting perturbation is the explicit sym- 
metry breaking between the dot-lead and dot-grain tun- 
neling amplitudes. To address this issue, it is convenient 
to rewrite the Kondo Hamiltonian in the most general 
form as follows (again r = for the bulk lead and r = 1 
for the grain): 



H K = 



E (JrriS-^ T 

T=0,1 ^ ' 
T=Q,1 ^ ' 



V 



+ -f[T + (^r-^) + h.c] 



(33) 



+ Q ± S- [T+{^T-ai))+h.c] . 
The corresponding bare values are embodied by 

Jo = J , Ji = K 2 J , Q± = 
V± = VK , Q Z . T = V Z . T = 0, 



(34) 



where we have introduced the asymmetry parameter 
K = ti/t; t = to (ti) denotes the hopping amplitude be- 
tween the lead (grain) and the small dot. Since the asym- 
metry stands for a marginal perturbation in the RG sense 



it is natural to argue that the SU(4) correlated ground 
state is still robust for weak asymmetry between the tun- 
nel junctions. But, to obtain more quantitative results we 
yet resort to NRG (Figure 10). By taking V± = 0.1 J, we 
can observe that the mixing of spin and orbital degrees 
of freedom may survive until K 0.95; This guarantees 
an anisotropy of roughly 10% between the conductances 
at the tunnel junctions to preserve the SU(4) fixed point. 
Mostly, the magnetic moment S and the isospin T are 
simultaneously quenched and again the spectrum can be 
understood as a sum of four independent chiral fermion 
spectra with phase shift 7r/4. 

Let us now discuss the case of a quite strong asymme- 
try between the tunnel junctions. For completeness, we 
also provide the RG equations at second order for this 
generalized situation 



dJ T 
IT 
dV z , T 
all 

dV± 

~dT 

dQ z , T 
dl 

dQ ± 



Jr 2 + {Q Z .r) 2 + 2Qj 



V± 2 + 3Q ± 2 



\v±_ (V z , + V x ,i) + (Q z , + Qz,i) 



dl 



2J r Q ZiT + 2V ± Q ± 



(35) 
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FIG. 10: Magnetic and orbital susceptibilities versus u>/D for 
close to unity values of the asymmetry parameter K — t/t\ 
between the two tunnel junctions. The SU(4) ground state 
is stable against the inclusion of a weak asymmetry between 
tunnel junctions. 

+ ^V± {Qz,a + Qz,i) ■ 

At second order, note the equality V Z ${1) — V Zt \{l) = 
V z (l) regardless of the parameter K. Primarily, it is 
immediate to observe that for K = 1, we recover the 
previous SU(4) Fermi-liquid flow. Now we greatly di- 
minish the tunneling amplitude between the grain and 
the small dot, i.e., t\ <C t (t being fixed) and K <C 1. 
With the present notations, it is clearly transparent that 
the spin Kondo coupling Jq = J between the bulk lead 
and the small dot will be the largest one through the 
RG flow and becomes of order unity at the temperature 
T K [K < 1] ~ De" 1/J = T K [A Z ] whereas all the other 
couplings are still negligible that breaks the SU(4) sym- 
metry explicity. 

It is worth noting at this stage that the role of the 




co/D 

FIG. 11: Spin susceptibility versus uj/D upon increasing the 
asymmetry between tunneling amplitudes at the two junc- 
tions. 

asymmetry parameter K seems to be practically equiv- 
alent to renormalize the orbital splitting A 2 (compare 
Figures 6 and 11). The main difference however is that 
at the degeneracy points of the grain, one can expect a 
second-stage quenching of the isospin T at some lower 
temperature, but obviously this (very) low-temperature 
regime lies much beyond the range of validity of the effec- 
tive Hamiltonian (J2HJ- Furthermore, one can clearly no- 
tice that the previous perturbative result of Eq. (|31(l di- 
verges if one of the charging energy U\OY U- i approaches 
zero, i.e., is not applicable. 

In fact, as already noted in Ref. JIH it is a very difficult 
task to find the exact shape of the step of the staircase 
in the present situation of a grain at a degeneracy point 
coupled to an effective resonant level. But qualitatively, 
one might expect that the physics and the resulting 
(two-channel) Kondo energy scale should not be so dif- 
ferent as those of a grain coupled to a normal lead with 
a reduced bandwidth Tj^[K <C 1], via a hopping matrix 
element t ~ t\: 

Tj? h =T K [K « lJe-T/' 1 ; (36) 

Here 7 is a constant parameter of the order of unity. 
A similar discussion should hold in the opposite regime 
K 3> 1 where one expect this time the Kondo coupling 
J\ to first flow to strong coupling (since it is proportional 
to K 2 ) at the temperature scale 7if[A 2 ] ~ De~ x / Jl . 
Since the conductance between the grain and the lead 
is still very small at the intermediate energy scale due 
to the anisotropy, a second stage quenching of the or- 
bital pseudo spin is expected at a lower energy scale in 
a similar manner as the case K <C 1- Unfortunately for 
asymmetric junctions, it is difficult to formulate more 
quantitative results at low temperatures. A complete 



13 



renormalization group calculation starting with the bare 
Hamiltonian would be necessary. This goes beyond 
the present analysis. Finally, let us mention that for more 
moderate values of U and e, i.e. rather in the resonant 
level regime of the small dot, the NRG results of Lebanon 
et alm^ still support a two-channel Kondo crossover and 
the ovserscreening of the isospin moment in the case of 
asymmetric junctions. 

D. Large junctions 

We predict that the SU(4) symmetry should be still 
robust for wider junctions characterized by n > 1 trans- 
verse channels with almost equal transmission ampli- 
tudes, however the associated Fermi-liquid typical energy 
scale decreases exponentially with the number of con- 
ducting modes. For instance, extending results of Ref. |2l| 
for our geometry, we can clearly assess that there will be 
a unique 'effective tunneling mode' in the lead (it is some 
combination of the original tunneling modes in the lead) 
and another unique 'tunneling mode in the box' (also a 
linear combination combination of the original modes in 
the grain): The T=0 effective Hamiltonian of the model 
at the degeneracy points of the metallic dot corresponds 
to tunneling between these two modes only with or with- 
out spin flip of the excess spin of the small dot, and all 
the other modes can be neglected. This entirely justi- 
fies the emergence of an SU(4) fixed point at very low 
temperatures even if the number of modes in the lead or 
in the grain is larger than one. However, the ultraviolet 
cutoff D at which the effective tunneling mode prevails, 
must be properly rescaled to^i 

T*[n] = De- an , (37) 

where a is of the order of unity. Unfortunately, this im- 
plies that an SU(4) Kondo singlet can only occur at the 
much reduced Kondo temperature scale 

Tf^^T'He- 1 '". (38) 

Experimentally, in order to maximize chances for observ- 
ing the SU(4) Fermi liquid realm, it is then more advan- 
tageous to consider tunneling junctions with one clearly 
dominant conducting transverse mode. 

VI. DISCUSSIONS AND CONCLUSIONS 

We have determined exactly the shape of the steps of 
the Coulomb staircase for a grain coupled to a bulk lead 
through a small quantum dot in the Kondo regime. First, 
we mapped the problem onto a related model of two ca- 
pacitively coupled small quantum dotsi2& Then, combin- 
ing both NRG calculations with perturbative scaling ap- 
proaches we have shed light on the possibility of a stable 
SU(4) Fermi liquid fixed point occurring at the degener- 
acy points of the grain, where a Kondo effect appears 
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FIG. 12: Sketch of the capacitance peaks for our setup with 
almost symmetric junctions (dashed line) compared to those 
in the original Matveev's problem (full line)i- 



simultancously both in the spin and the orbital sectors; 
This demands symmetric or slightly asymmetric tun- 
nel junctions and preferably a single-conducting chan- 
nel with two spin polarizations. More generally, as in 
Ref. pr|. these results bring precursory insight on the 
realization of Kondo ground states with SU(N) (N = 4) 
symmetry at the mesoscopic scale. 

Let us provide a physical interpretation for the occur- 
rence of such an SU(4) entanglement. Typically, close 
to the degeneracy points of the grain, we have two spin 
objects, namely the spin S of the small dot and the or- 
bital pseudo-spin T of the grain depicting the two al- 
lowed degenerate charging states. Obviously, when these 
two spin objects are uncoupled the symmetry group of 
the problem is unambiguously SU{2) ® SU(2). But, as 
already discussed at length in the sequel, in our setting 
spin- flip assisted tunneling events - i.e., an electron from 
the bulk lead tunnels onto the metallic grain by flipping 
the excess spin of the small dot and vice- versa - are very 
prominent at low energy; This implies that the infra-red 
fixed point must also reflect a visible spin-orbital mix- 
ing. Finally, it is easy to check that SU(4) is the mini- 
mal group allowing spin- orbital entanglement and which 
guarantees rotational invariance both in spin and orbital 
spaces. Our Kondo fixed point then is rather described 
by the quenching of the hyper-spin [S + |][T + i]. 

In a very different context, let us mention that SU(4) 
singlets have also shown up in fermion lattice models 
where spin and orbital degrees of freedom play a very 
symmetric rolei^SiiSi 

The major consequence of this enlarged symmetry is 
that the ground state is Fermi-liquid like, which consider- 
ably smears out the Coulomb staircase behavior already 
in the weak tunneling region, and in particular, hinders 
the Matveev logarithmic singularity 9 to take place (Fig- 
ure 12). The grain capacitance exhibits instead of a loga- 
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rithmic singularity, a strongly reduced peak as a function 
of the back-gate voltage. This stands for an irrefutable 
signature of the formation of a Fermi-liquid ground state 
when tunneling through a single- electron box. Further- 
more, we mightily emphasize that our NRG calculations 
markedly reproduce an SU(4) ground state regardless 
of the particle-hole asymmetry onto the small dot (Fig- 
ure 4); More precisely, even in the case of particle-hole 
symmetry 2e + U = 0, the spectrum can be still inter- 
preted as a sum of four independent chiral fcrmions with 
phase shift 7r/4 in agreement with the SU(4) Fermi liq- 
uid theory. This differs from the conclusion of Ref . |23 • 
However, this is not so surprising in the sense that in 
their NRG calculations (see, e.g., their Figures 15 and 
16), Lebanon et al. have studied a rather different limit, 
U = — 2e but U/E c <C 1, which does not correspond to 
our situation of a small dot and a much larger metallic 
grain (U/E c ^S> 1). Besides, in the case of symmetric 
barriers, they clearly noticed that a moderate Coulomb 
repulsion on the small dot already pushes the two-channel 
Kondo regime down to much lower temperature. 

It is also worth to recall that the associated Kondo 
temperature scale T^ U ^ can be strongly enhanced com- 
pared to that of the Matveev's original setup which 
maybe ensures the verification of our predictions. In 
particular, for very large U [U — 2e and V± > 0], 
exp— (1/4 J) may be larger than the Kondo 
scale in the conductance experiments across a single small 
quantum dol£ (~ IK), and capacitance measurements 
can be performed much below lOOmKmll Additionally, 
we have checked that the SU(4) Kondo temperature scale 
is considerably reduced for negative values of V±, i.e., 
upon by (moderately) decreasing the on-site interaction 
U (U ~ — e), i.e., by making the small dot larger and 
larger4i We have carefully discussed the robustness of 
the SU(4) correlated state against the inclusion of weak 
perturbations like an external magnetic field, a deviation 
from the degeneracy points, or still an asymmetry in the 
tunnel junctions. 

Let us now pursue and discuss an interesting crossover. 
So far, we have concentrated on the situation at and near 
the degeneracy points of the grain. Let us now apply a 
quite strong orbital magnetic field such that we move ex- 
plicitly away from the degeneracy points. Naively, since 
one suppresses the orbital spin-flip terms, one could in- 
fer the emergence of a two-channel spin Kondo model 
through the two Kondo terms Jo and J\\ However, in our 
setting with almost symmetric junctions, the Schrieffer- 
Wolff transformation away from the degeneracy points 
always ensures Jo > Ji; The NRG calculation of Fig- 
ure 6 clearly reproduces this expectation. The system 
then undergoes a one-channel Kondo crossover. First, 
the emergence of a logarithmic contribution in (Q) at 
quite high temperature could be potentially observable. 
Furthermore, at low energy, the physics resembles that 
of a resonant level - induced by the formation of an 
Abrikosov-Suhl resonance between the small dot and the 
bulk lead - weakly-coupled to the grain; We then recover 
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FIG. 13: Another mesoscopic double lead setup, candidate 
for the SU(4) model. This could be equally performed with 
vertically coupled dotsi— 



a similar situation to that of Ref. [24j . 

Another possible realization of our SU(4) model could 
be still possible in a multi-lead geometry (Figure 13). 
Again, this would demand to be at the degeneracy points 
of the grain and to adjust the different tunneling junc- 
tions. More precisely, following Glazman and Raikh, 34 
only the even linear combination of the electron creation 
and annihilation operators in the two bulk leads couples 
to the local site (small dot). The odd linear combina- 
tion can be omitted and conceptually the effective model 
could be rewritten as in Eq. I|17fl . Let e.g. assume that 
the tunnel junctions between each lead and the small 
dot are symmetric. Then, only the linear combination 
"00 = (V'oi + ip02)/V^ will be coupled to the small dot; 
tpoi (i = 1,2) denotes the electron annihilation operator 
in each lead. To recover an SU(4) Kondo fixed point, we 
infer that the grain-dot tunneling amplitude then must 
be approximately V2 times that between each lead and 
the small dot. This setup is particularly interesting be- 
cause the capacitance of the grain and the conductance 
across the small dot could be both measured. Further- 
more, blocking completely the opening between the grain 
and the small dot, one could recover a more usual Fermi 
liquid behavior with SU(2) spin symmetry when measur- 
ing the conductance across the small dot, and observe a 
net reduction of the Kondo energy scale compared to the 
SU(4) case due to spin orbital decoupling. 

Note that this geometry - away from the degeneracy 
points of the grain - has been previously discussed by 
Oreg and Goldhaber-Gordon as a potential candidate for 
the appearance of a two-channel (spin) Kondo regime in 
a conductance measurement^ This requires meticulous 
fine-tuning of the gate voltages and tunnel junctions to 
equalize the coupling to the two channels (grain plus even 
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linear combination of the leads) . 

Definitely, the potential observation of a two-channel 
Kondo effect in artificial nanostructures would be an im- 
portant issue2^iS4i since the emergent non-Fermi liq- 
uid behavior is very intriguing and so far difficult to ob- 
serve with real magnetic impurities due to the intrinsic 
channel anisotropy. 15 In our setting, another interesting 



experiment to do in order to have potentially access to 
a two-channel (charge) Kondo behavior would be to stay 
at the degeneracy points of the grain and then progres- 
sively to shift the impurity level e on the dot (which can 
be tuned via the gate voltage Vd of the small dot) to the 
Fermi energy, i.e., to reach the mixed- valence (=resonant 
level) limit for the small dot~ 
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APPENDIX A: PERTURB ATIVE CALCULATIONS 



Here, we derive explicitly the perturbative result of Eqs. (JSJ and (9). We essentially focus on the Kondo term; the 
perturbation theory for the direct hopping term V can be found in Ref. First, it is accurate to rewrite the Kondo 
term in real space as: 



H, 



E 

a0 



Ju 



3=0,1 



(Al) 



where ipo a = J2k a ka and "01a = J2 P a pa- The granule charge operator reads Q = eJ2 a V'laV'ia- Now, let |0) denote 
the ground state of the unperturbed Hamiltonian with t — — oo. The first order correction |1) to |0) then readst^ 



11) = -* 



dtH K {t)\0), 



(A2) 



Hk being taken in the interaction representation. The expectation value of the charge on the dot however is second 
order in the Kondo coupling; Indeed, we easily get (0|Q|1) = 0. Therefore, the most leading contribution takes the 
form {Q)i = (0|Q( 1 "'|1), where is the first order correction to the charge operator on the dot. This can be 
computed using the identification: 



J(t)dt with J(t) = i[H K , Q] 



J must be identified as the effective current operator mediated by the Kondo coupling. This results in 



Q 



(i) _ 



J 



a/3 



dt 



The expectation value of the charge on the dot is then to second order in the coupling to the impurity 



(Qh 



(Ji. 



EE 



a, b aft 



dt% 



3(Jio) 2 



dt x 



dU 



dt 2 (5 a (ti)^(t 2 )aV b )[(^ Q (^)^oa(ii))(V'i«(i2)V'L(ii)> 

X) 



(A3) 



(A4) 



(A5) 



where the averages are taken over the ground state of the uncoupled system. It is advantageous to Fourier transform 
the problem as: 



(Q>2 = - 



3(Jio) 2 



E 

p,k 



dti 



dU 



(afe(i 2 )4(*i))( a p( < 2) a p( i i)> ~ M*2)4(ti))(4(*2H(ti)) 



(A6) 
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where the momentum indices p and k respectively refer to the grain and to the reservoir. Using the Green's functions 
of the isolated grain: 

(al(t 2 )a p ( tl )) = eC-epy^-^Jte-'O (A7) 
(a p (t 2 )al( tl )) = e(e p )e-^ +u ^-^, 

where again U\ and U—i embody the energies to add an electron and hole onto the grain, we finally find: 

3(Jio) 2 W e(6 fc )6(-£ p ) 6(- efc )6(e p ) 1 g 



(Q)2 = -e — - — 2^ 



(e fe -e p + U-xY (e p - e k + f7i) 2 



3(Jio) 2 , (e/2C-ip 
m 1 



8 \e/2C + <p 

is the usual Heavyside function. Density of states in the grain and in the lead have been assumed to be equal and 
taken to be 1 for simplicity. 

Now, we briefly want to show that cubic orders involve logarithmic divergences both associated with the Kondo 
coupling and with the proximity of a degeneracy point in the charge sector. More precisely, let us focus on the specific 
contribution in Jo</ 2 for the term (Q) 3 = (0\Q^\2), with 

1 f° r° 

|2} = -- / dt x / dt 2 T[H K (ti)H K (t 2 )]\0) (A9) 

* J-oo J — oo 



JqJ: 
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EEE f *i f d^nS%h)S\t 2 )}T[4 a {t 1 )^off(tx)^(t 2 )^ iM*a)]|0) 

a,b a/3 n,v J ~°° J ~°° 

EEE f *i f ^2r[^(ti)^fe)][r(^L(ii)^(t 2 ))^,^(ii)v4^)^^]|o) 

a, 6 a/3 n,u J -°° 

EEE/ / * 2 5 G S5 n(t 1 -t 2 )T(vL(ii)V'o a (t2)}Vl M (*2)^Vo/3(ti)|0) 



2 

c a jj, t f) 



J oJioEE ln (VV) / ^U(*i)^^)/j(*i)|o). 



c fj.,/3 

It becomes then obvious that |2) is (almost) proportional to |1); It is straightforward to show that this induces a 
third-order correction for the charge on the grain 

(Q(T)>3OcJ (J 10 fln(^)ln(^ T |). (A10) 

Note that the appearance of the extra In (D/T) factor clearly stems from the prominent renormalization of the lead-dot 
spin Kondo coupling Jq on a charge plateau. 
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